Nernst effect beyond the relaxation-time approximation 
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Motivated by recent interest in the Nernst effect in cuprate superconductors, we calculate this 
magneto-thermo-electric effect for an arbitrary (anisotropic) quasiparticle dispersion relation and 
elastic scattering rate. The exact solution of the linearized Boltzmann equation is compared with 
the commonly used relaxation-time approximation. We find qualitative deficiencies of this approxi- 
mation, to the extent that it can get the sign wrong of the Nernst coefficient. Ziman's improvement 
of the relaxation-time approximation, which becomes exact when the Fermi surface is isotropic, also 
cannot capture the combined effects of anisotropy in dispersion and scattering. 

PACS numbers: 72.15.Jf, 73.50.Jt, 74.25. fc, 74.72. Gh 



I. INTRODUCTION 

The Nernst effect is a magneto-thermo-electric ef- 
fect, in which an electric field in the a;-direction 
results from a temperature gradient dT/dy in the y- 
direction, in the presence of a (weak) magnetic field 
B in the z-directionii The Nernst coefficient Mxy ~ 
—Ea:{BdT/dy)~^ depends sensitively on anisotropics in 
the band structure. In particular, for a square lattice 



xy 



■Afyx is antisymmetric upon interchange of x and 
y — just like the Hall resistivity — but lattice distortion 
breaks this antisymmetry. 

There has been much recent interest in the Nernst ef- 
fect in the context of high- Tc superconductivity, since 
underdoped cupratcs were found to have an unusually 
large Nernst coefficient in the normal state.— This may 
be due to superconducting fluctuations above Tj,, or it 
may be purely a quasiparticle effect The quasipar- 
ticle Nernst effect has been studied on the basis of 
the linearized Boltzmann equation in the relaxation- 
time approximation.—"— This is a reliable approach 
if the scattering rate is isotropic, since then the ne- 
glected "scattering-in" contributions average out to zero. 
There is, however, considerable experimental evidence 
for predominantly small-angle elastic scattering in the 
cupratesr^"— possibly due to long-range potential fluc- 
tuations from dopant atoms in between the Cu02 
planes.—"^ 

It is not surprising that existing studies rely on the 
relaxation-time approximation, since the full solution 
of the Boltzmann equation with both band and scat- 
tering anisotropics is a notoriously difficult problem.— 
In our literature search we have found magneto- 
electric calculations that go beyond the relaxation-time 
approximation^i^"— but no magneto-thermo-electric 
studies. It is the purpose of this paper to provide such a 
calculation and to assess the reliability of the relaxation- 
time approximation. 

We start in Sec. |ll] with a formulation of the anisotropic 
transport problem, in terms of the socalled vector mean 
free path.— In the relaxation-time approximation, this 
vector Afc is simply given by the product VkTk of veloc- 
ity and scattering time (all quantities dependent on the 



point fc on the Fermi surface). Going beyond this ap- 
proximation, Afc is determined by an integral equation, 
which we solve numerically. 

We also consider, in Sec. IIIIl an improvement on the 
relaxation-time approximation, due to Ziman,— which 
incorporates some of the scattering-in contributions into 
the definition of the scattering time. For isotropic Fermi 
surfaces Ziman's scattering time is just the familiar trans- 
port mean free time — which fully accounts for scatter- 
ing anisotropics. If the dispersion relation is not isotropic 
this is no longer the case. 

Wc compare the exact and approximate solutions in 
Sec. llVland conclude in Sec. IVl 



II. FORMULATION OF THE TRANSPORT 
PROBLEM 

A. Boltzmann equation 

We start from the semiclassical Boltzmann transport 
equation for quasiparticles (charge e) in a weak magnetic 
field B, driven out of equilibrium by a spatially uniform 
electric field E and temperature gradient VT. The exci- 
tation energy is e^, relative to the Fermi energy ep- The 
band structure may be anisotropic, so that the velocity 



Vk 



(2.1) 

(with Vfc = d/dk) need not be parallel to the momentum 
hk. For simplicity, we assume there is only a single type 
of carriers at the Fermi level (cither electrons or holes). 

Upon linearization of the distribution function fk = 
fo + 5fc around the equilibrium solution 

Jo- 



l + cxp[{ek-eF)/kBTy 
the Boltzmann equation takes the forroii 
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E 



Q{k,k'){gk 



U 



eE 



£fc — £f 
T 



(2.2) 

9k'), 
(2.3) 
(2.4) 



2 



The right-hand-side of Eq. p.3|) is the difference be- 
tween the scattering-in term Q{k, k')gf,i and the 
scattering-out term J2k' Qi^' ^^)9k (with Q{k\k) ~ 
Q{k, k') because of detailed balance). 
We assume elastic scattering with rate 

Q{k,k')^d{ek~ek')q{k,k') (2.5) 

from fc' to fc. Detailed balance requires 

q{k\k)^q{k,k') (2.6) 

and particle conservation requires 

E-9fe = 0- (2-7) 
fc 

The sum over k represents a d-dimensional momentum 
integral, J2k ^ (Stt)"''/ dk ( in a unit volume) . The spin 
degree of freedom is omitted. 

It is convenient to define the Fermi surface average 



■dSFf{k)\Vk\ 
USF\Vk\-' 



(2.8) 



with a weight factor \vk\ ^ from the volume element dk 
fi^^lvkl^^dekdSp- The density of states is given by 

N {e f) ^ h- \2Tr)~ <l^dSF\vk\'\ (2. 

For later use we note the identity 



9) 



{fik)ivk X Vfc)5(fe))5, - -{9ik){vk X Vk)f{k))s,, 

(2.10) 

valid for arbitrary functions /, g of fc. 



They can be written in terms of Fermi surface averages, 
N{eF){q{k, fc')(Afc - Afe'))s^ = Vk, (2.14) 

N{eF){q{k,k'){6Ak - SAk'))s' = ■ {vu x Vk)^k- 

(2.15) 

(The prime in the subscript S'p indicates that k' is av- 
eraged over the Fermi surface, at fixed fc.) The solution 
should satisfy the normalization 



(Afc)s, = = {8Ak)s, 



(2.16) 



required by particle conservation to each order in B. 

The integral equations (|2.12p and (|2.13p can be readily 
solved numerically. In the limit of small-angle scattering 
an analytical solution is possible, by expanding the fc'- 
dependence around fc to second order,— but we have 
not pursued that method here. 



C. Linear response coefficients 

In linear response the electric current density j is re- 
lated to the electric field E and temperature gradient VT 

by 



j =crE aVT. 



(2.17) 



The conductivity tensor cr follows from the vector mean 
free paths by 



dgk 
BE 

dfo 



J2[-g^)vk®i^k + 5Ak) 



(2.18) 



B. Vector mean free paths 



We seek the solution of Eq. (|2.3p to first order in B. 
Following Refs. [22II23I we introduce the vector mean free 
paths Afe (of order 5°) and 6Ak (of order B^), by sub- 
stituting 



gk = U ■ (Afc + SAk) . 



(2.11) 



Since the vector U can have an arbitrary direction it 
cancels from the equation for Afc. The equation for 6Ak 
has also a term oc {vk x Vk)U, which vanishes because 
Vfci7 = hvkdU/dek- 

The resulting equations for the vector mean free paths 
are 



Q{k, fc')(Afc - Ak') = Vfc, 



(2.12) 



^ Q(fc, fe')(^Afe - 5Ak') = |b • {vk X Vfc)Afc. 



(2.13) 



[The direct product indicates a dyadic tensor with ele- 
ments (a (g) b)ij = ttibj.] 

At low temperatures, when —dfo/dsk — !> S{ek — Sf), 
this may also be written as a Fermi surface average. 



cr = e^N{eF){vk ® (Ak + S Ak)) s ■ 



(2.19) 



By substituting Eq. p.l4|) for Vk and using Eq. (|2.15|) 
together with the detailed balance condition (|2.6p and 
the identity (|2.10p . one verifies the Onsager reciprocity 
relation 

a,,(B) =a,,(-B). (2.20) 
The thermoelectric tensor a is given by 
dgk 



a 



fc 

e 
f 



d{-\7T) 



dfo 



^iSk - £f) (^^^ ) ® (^'^ + '^^fc) • 



(2.21) 
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At low temperatures this reduces to the Mott formula, 

(2.22) 



_B-independent Nernst coefficients 



a = ^ — 



3e dep 

These equations all refer to a single type of carriers 
at the Fermi level (electrons or holes), as would be ap- 
propriate for hole-doped cuprates. The ambipolar effects 
of coexisting electron and hole bands are not considered 
here. 



D. Nernst effect 

We take a two-dimensional {d = 2) layered geometry 
in the x — y plane, with a magnetic field B = Bz in the z- 
direction. The Nernst effect relates a transverse electric 
field, say in the a;-direction, to a longitudinal temperature 
gradient (in the y-direction) , for zero electric current. 

One distinguishes the isothermal and adiabatic Nernst 
effectji depending on whether dT/dx = or jh,x = is 
enforced (with jh the heat current). As is appropriate 
for the cuprateS)^ we assume that a high phonon contri- 
bution to the thermal conductivity keeps the transverse 
temperature gradient dT/dx negligibly small, so that the 
Nernst effect is measured under isothermal conditions. 

The isothermal Nernst effect is expressed by 



xy/B 
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Myx = - 



3e 



^ hm 



B^b B Gxx dep cTyy 

1 (Txx d, (T 
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3e B^o B <jyy dep a. 



(2.26a) 
(2.26b) 



These expressions relate the Nernst coefficients to the en- 
ergy derivative of the Hall angle in the small magnetic- 
field limit. The cancellation in Eq. (j2.26ap of any iden- 
tical energy dependence of Gxy and Uyy is known as the 
Sondheimer cancellation»^i2^ On a square lattice one has 
(^xx = CTyy, hence Mxy = —Myx^ but without this C4 
symmetry the two Nernst coefficients differ in absolute 
value. 

In terms of the vector mean free paths, the Nernst 
coefficients arc given by 



^ _ TT kgT (vk,y^k,y}sF ^ {vk,x5Kk,y) Sf 

'ieB (vk,x^k,x)sF dep {vk,yA.k,y)sF 



(2.27a) 



^ _ 7r^fc|r {vk,x^k,x)sF d {vk.xSAk,y)sF 
3eB {vk.yAk,y)sF dep {vk,xAk,x)sF 

(2.27b) 

where we have used that A;, is _B-independent and 6Ak 
is oc B. 



Ex — Ox 



dT 
dy 



dT 
dx 



0, 



and similarly with x and y interchanged, 
mopower tensor 



has off-diagonal elements 



(2.23) 



The ther- 



(2.24) 



III. RELAXATION-TIME APPROXIMATION 

In the relaxation-time approximation the scattering-in 
term Q(fe, fc').9fc' the right-hand-side of the Boltz- 
mann equation p.3|) is omitted.— Only the scattering-out 
term gk 'Yliw QC'^j '"O ~ 9k/ Tk is retained, containing the 
momentum dependent relaxation rate 



l/Tfe = ^Q(fc,fc') = N{ep){q{k,k'))s'^ 



(3.1) 



Jyx 





O'xyO^yy 


O'xx^yy 
^xx^yx 


O'xyO'yx 

O'yxOtxx 


^xx^yy 


— (JxyO'yx 



(2.25a) 
(2.25b) 



We will consider two-dimensional anisotropic band 
structures that still possess at least one axis of reflection 
symmetry, say the y-axis. Upon reflection the component 
I— > —jx of the electric current changes sign, while Ey 
and dT /dy remain unchanged. The perpendicular mag- 
netic field B ^ —B also changes sign, because it is an 
axial vector. It follows that axy{B) = —axy{—B) and 
ctxy{B) = —axy{—B) are both odd functions of B, so 
they vanish when _B — > 0. 

Using the Mott formula (j2.22p , one can then define the 



Without the scattcring-in term, the equations (j2.12p 
and (|2.13p for the vector mean free paths can be solved 
immediately, 

e 

Afc = DfcTfc, 5Ak = J^'^kB ■ {vk X Vfe)Tfct>fc. (3.2) 

In general this solution does not satisfy the particle con- 
servation requirement (|2.16|) . which is the fundamental 
deficiency of the relaxation-time approximation. 

Substitution into Eq. ([239]) gives the conductivity ten- 
sor 



cr = e^N{ep){TkVk ® {vk + ^kTkVk))sF^ 
with differential operator 



-B ■ {vk X Vfc) 



(3.3) 



(3.4) 
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For a two-dimensional lattice with reflection symmetry 
in the y-axis, the elements of the conductivity tensor are 
given by 

cTxx = e^N{eF){Tvl)sp, (Tyy = N {e f) {tvI) s p , (3.5) 



d d 

X ( TV,,, I Vxjrr '^'v^rT ] '^'"y 

UKy OKx / / Sf 



(3.6) 



(Here we don't write the subscript k to simplify the no- 
tation.) The Nernst coefficients in the relaxation-time 
approximation then follow from Eq. (|2.26|) as the energy 
derivative of the ratio of two Fermi surface averages, 



^yy 



J^xy — Zq- 

a XX «£/ 



TVx [vx 



at „ '"y at^ } ''""y 



(3.7a) 



yx 



-z„ 



TVx {vx 



Gyy 



where we have defined 



(3.7b) 



(3.8) 



One may further simplify the relaxation-time approx- 
imation by taking an isotropic relaxation time tq^ef), 
which is the approach taken in Refs. IcI-ITol Since {vk x 
Vfc)To(£F) = 0, Eq. p.7p then reduces to 



J^xy — ^0 



CTxx dSF (v: 



y/SF 



, dv,, 



dv„ 



VxVy 



cfxx d ro(£F) I ^,2^'^y ^"y 



. dv„ 



dh 

dv. 



(Jyy dSF {vDSf X""" dky ^""^^ Otx 



(3.9a) 



Sf 

(3.9b) 



If one stays with a momentum dependent relaxation 
time Tfc , then it is possible to improve on the relaxation- 
time approximation by changing the definition (j3.ip into 
Ziman's expressio n^^i^^ 



l/rf--=iV(£F^)(g(fe,fc')(l- 



Vk ■ Vk' 

\vk\ \vk'\y I 5^ 

('3.IO) 

Ziman's improvement of the relaxation-time approxima- 
tion becomes exact if the Fermi surface is isotropic, mean- 
ing that Ek is only a function of |fc| and q{k,k') is only 
a function of fc • fc'. 



IV. COMPARISON 

We turn to a comparison of the Nernst effect in 
relaxation-time approximation with the exact solution 
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FIG. 1: Dependence of the Nernst coefficient on the distortion 
e of the square lattice at a fixed band filling riband = 0.875, for 
three different values of the range ^ of the scattering poten- 
tial. The three panels show how the exact solution of the lin- 
earized Boltzmann equation (solid) starts out very close to the 
relaxation-time approximation (dotted) for nearly isotropic 
scattering, and then becomes progressively different as small- 
angle scattering begins to dominate. 



of the linearized Boltzmann equation. For this com- 
parison we need to specify an elastic scattering rate 
Q{k,k') = 5{£k — £k')q{k,k') and a dispersion relation 
£k- 

For the scattering, we take a random impurity poten- 
tial with range ^. By increasing ^ relative to the Fermi 
wave length, we can study the transition from isotropic 
scattering to (small-angle) forward scattering. We model 
the impurity potential by a sum of Gaussians, centered 
at the random positions r.^ of the impurities, 



U{r) 



Ui cxp 



\r - n| 

e 



(4.1) 



The amplitude Ui is uniformly distributed in [—(5, (5]. The 
correlator is 



(C/(r)f/(r'))--<5TriimpCxp 



,.'|2 



2e 



{\U{k)\ 



l^^^^Vmp exp(-iC2|fc| 



(4.2) 
(4.3) 



where nimp is the two-dimensional impurity density 
(number of impurities per area per layer) . The resulting 
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FIG. 2: Dependence of the Nernst coefficient on the range 
^ of the scattering potential, for an undistorted square lat- 
tice (e = 0). Two values of the band filling are shown in the 
upper and lower panel. The three curves in each panel cor- 
respond to: the exact solution of the linearized Boltzmann 
equation (solid), the relaxation-time approximation (dotted), 
and Ziman's improvement on the relaxation-time approxima- 
tion (dash-dotted). 



elastic scattering rate (in Born approximation) becomes 
q{k,k') =7oexp(~ie'|fc-fc?), 

70 = — -■ (4.4) 



6h 

Values of ^/a of order unity are to be expected in the 
cuprates for scattering by impurities between the Cu02 
planes, when ^ is of the order of the interplane distance. 

For the dispersion relation we follow a recent study 
of the Nernst effect in hole-doped cupratesji by taking 
the tight-binding dispersion of a distorted square lattice 
with first (ti), second (^2), and third (^3) nearest-neigbor 
hopping: 

E{k) = - 2ti [(1 + e) cos fc^; + (1 - e) cos ky] 
- 2*3 [(1 + e) cos 2k^ + (1 - e) cos 2fc„ 



4*2 cos kx cos ky. 



(4.5) 



The lattice constant is a and k is measured in units of 
1/a. The C4 symmetry is distorted by the anisotropy 
parameter e, preserving reflection symmetry in the x and 
y-axes. 
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FIG. 3: Same as Fig. [2] but now showing the dependence on 
the distortion e of the square lattice for a fixed range ^ = 
0.75 a of the scattering potential. 



We use ratios of hopping parameters t2/*i = 0.32, 
^3/^2 = 0.5, and compare two values of the band filling 
fractions riband = 1.156 and 0.875. (The corresponding 
Fermi energies at e = are Ep = and Ep w — 0.97*i 
respectively, and are adjusted as e is varied to keep riband 
fixed.) 

The Nernst coefficient is plotted in units of 



A/'o 



(4.6) 



'iK^ 70 

We show only JVxy, since Myx is related by 

Uxy{e) = -Myxi-e). (4.7) 
We compare three results for the Nernst coefficient: 

• the exact solution of the linearized Boltzmann 
equation, from Eq. p.27|) : 

• the momentum-dependent relaxation-time approx- 
imation, from Eq. p.7p : 

• Ziman's improvement on the relaxation-time ap- 
proximation, from Eq. p.lOj) . 

We have have found that there is little difference between 
the momentum-dependent and momentum-independent 
relaxation-time approximations [Eqs. p.7p and p.9p ]. so 
we only plot the former. 
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Results are shown in Figs. [T]-[3l Fig. [T] shows that the 
relaxation-time approximation agrees well with the exact 
solution for nearly isotropic scattering <C a). With 
increasing ^ small-angle scattering begins to dominate, 
and the relaxation-time approximation breaks down for 
^ > 0.4 a. In Fig.[2]we see that Ziman's improved approx- 
imation remains reliable over a somewhat larger range of 
^. Still, for a modestly large f = 0.75 a also Ziman's 
approximation has broken down completely, see Fig. |3j 
giving wrong magnitude and sign of the Nernst coeffi- 
cient. 



V. CONCLUSION 

In conclusion, we have shown that the relaxation-time 
approximation is not a reliable method to calculate the 
Nernst effect in the combined presence of band and scat- 
tering anisotropics. The deficiencies are qualitative, even 
the sign of the effect can come out wrong. Of course, the 



relaxation-time approximation remains a valuable tool 
to assess the effects of band anisotropy in the case of 
isotropic scattering. 

Wc have based our comparison on parameters relevant 
for the cupratesji but we have only considered one pos- 
sible mechanism (single-band elastic quasiparticle scat- 
tering) for the Nernst effect in cuprate superconductors. 
Other mechanisms (ambipolar diffusion, inelastic scat- 
tering, superconducting fluctuations) would require sep- 
arate investigations!^ It is hoped that the general frame- 
work provided here will motivate and facilitate work in 
that direction. 
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